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The first-order interaction correction to the irreducible polarization function of pristine graphene 
is studied at arbitrary relation between momentum and frequency. The results are used to cal- 
culate the dielectric function and the dynamical conductivity of graphene beyond the standard 
random-phase approximation. The computed static dielectric constant compares favorably with 
recent experiments. 
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I. INTRODUCTION 

The influence of Coulomb interactions on elect ron 
properties of graphene is a subject of much interest 
Prominent examples of interaction effects are the en- 
hancement of the quasiparticle velocity v above the typi- 
cally quoted value v = 1.0 X 10 8 cm/s at low momenta^HSl 
(as predicted by the early theoretical worle^h anc j the 
emergence of additional dispersion branches^ believed 
to be plasmarons — bound states of electrons and plas- 
mons. The strength of Coulomb interactions in graphene 
is determined by the dimensionless parameter (h = 1 in 
this paper) 
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which can be controlled by varying the effective dielectric 
constant n of the graphene environment. 

Interactions modify not only the quasiparticle proper- 
ties but also the response functions of graphene, e.g., 
the irreducible (or proper) polarization P(q,uj). This 
quantitiy is defined to be the coefficient of proportion- 
ality between the change in electron density and the self- 
consistently determined screened scalar potential. The 
observables directly related to P are the dielectric func- 
tion e and the longitudinal conductivity a: 
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e(«,wj = 1 P{q,u), 
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a{q,u) = e -j P(q,uj) 
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Our work is motivated in part by two recent 
experiment j 10 ! 11 ! that suggested that the standard ran- 
dom phase approximation 12 (RPA) significantly underes- 
timates the static dielectric function e(q, 0) of graphene. 
The RPA amounts to replacing the exact polarization 
function P by the noninteracting value Pq. For neutral 
graphene at zero temperature, which is the case studied 
here, function Pq is given bj0 
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FIG. 1. First-order diagrams for the polarization function in- 
clude the self-energy part (a) and the vertex part (b). Solid 
lines represent electron Green's function, dashed lines the 
Coulomb interaction, and dotted lines the external momen- 
tum and frequency. 



at small enough q and oj where the Dirac approximation 
is valid. Equations ^ and Q entaiP^ 



eRPA(<7,0) = 1 + - a. 
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For k = 1 this formula gives crpa ~ 4.6. On the other 
hand, a much larger value e = 15.4^g 9 4 6 was inferred 
from inelastic x-ray scattering on bulk graphiteP^ Simi- 
larly, for graphene on a boron nitride substrate (k 2.5, 
a ps 0.9) crpa ~ 2.4, whereas the study of charge profile 
near Coulomb impurities by means of scanning tunneling 
microscopy suggests^ e = 3.0 ± 0.1. 

The first-order interaction correction P\ = O(o) be- 
yond the RPA is represented by the diagrams depicted 
in Fig. [T] We show that by including this correction, 



P(q,u) -> P a (q,cu) + P 1 {q,uj) 
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one can significantly reduce the discrepancy between the 
theory and experiment. Our result for e(q, 0) ii^ 



e(q,0) = 1 + — a q 
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(The difference between a q and a will be clarified in 
Sec. |nj) This formula yields e « 3.0 at a q = 0.9, in ex- 
cellent agreement with Ref. [TTJ In turn, if we use Eq. 
for a q — 2.2, we get e f=s 8.2, which is close to the lower 
estimate of the static dielectric constant in Ref. [10l al- 
though, the use of a perturbative formula is questionable 
at such large interaction strengths. 
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Another topic of interest is the dynamic response of 
graphene. For example, in the "optical" limit oj vq, 
Eqs. ([3]) and Q imply that noninteracting Dirac fermions 
have a frequency-independent conductivity^ 
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Accordingly, deviations of the infrared optical conduc- 
tivity of graphene from the universal value 00 may signal 
interaction effects. Exper imentally, no such deviations 
have been observe d 4 * 16 * 17 ! while the theoretical calcula- 
tion of the interaction corrections has been a subject of 
a debate. Our analysis favors the result 
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derived in Refs. [TS] and We explain why a different 
coefficient was obtained in Ref. [301 The numerical small- 
ness of (19 — 67r)/12 ps 0.01 may be the reason why the 
interaction corrections have not been observed so far. 

To get a more complete understanding of the inter- 
action corrections we also compute Pi(q, 0;) at arbi- 
trary Lu/vq ratios. This enables us to see the evolu- 
tion of Pi(q,u)) in the full momentum-frequency param- 
eter space, from the small value in the optical limit 
to the divergence at the spectral boundary uj — vq to 
a sizable effect in the static limit. Besides theoretical 
interest, motivation for this calculation comes from a 
rapidly burgeoning effort in probing regimes of interme- 
diate u>/vq by state-of-the-art experimental techniques, 
such as the near-field optics^M] and electron energy loss 
spectroscopies!^!!] Such regimes are also pertinent for 
the electromagnetic response of graphene nanoribbonsP^ 

The remainder of the paper is organized as follows. 
In Sec. HI] we summarize our main results. Derivation of 
these results is outlined in Sec. |III| In Sec. |IV| we compare 
our findings with previous work. The Appendix [X] is de- 
voted to the critical analysis of the controversy regarding 
the optical limit. 



II. MAIN RESULTS 

In order to present our results we first need to explain 
the notation a q in Eq. Q above. Crudely, v and a rep- 
resent the bare velocity and the bare coupling constant 
of the theory, whereas v q and a q denote their renormal- 
ized values. More precisely, v q is defined to be the phase 
velocity of quasiparticles with momentum q: 
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where S g is the on-shell self-energy. At the level of the 
first-order perturbation theory one 
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where A is the high-momentum cutoff. Therefore, v q and 
a q are given by 
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According to the renormalization group approach, such 
expressions are not just first-order approximations. They 
are asymptotically exact at low enough q where a q <C 1. 
However, a should be understood as the running coupling 
constant evaluated at some other cutoff A. The choice of 



a is largely arbitrary because a change in a in Eq. ( 13 ) 



can be absorbed into A. Nevertheless, a q , which is deter- 
mined by the observable quantity, the phase velocity v q , 
is unambiguous. Relations between various observable 
quantities are expressible in terms of the renormalized 
parameters only. For example, phase velocities at two 
different momenta q and k are linked by the relation 
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which is free from the arbitrary parameters a and A. 
Similarly, the polarization function can and (if one desires 
higher accuracy) should be expressed in terms of a q and 
V q . This point will be stressed again in Sec. IIVI 



Let us now present our findings. The first-order cor- 
rection to the polarization function is written as 



Pi(q,u) = —pi(x) 

Vq 



where x is the dimensionless ratio 
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and pi(x) is the complex dimensionless function, whose 
real and imaginary parts are displayed in Figs.[2ja)-|2|d). 
In general, pi[x) has to be evaluated numerically. How- 
ever, analytical results are available in several limits, as 
discussed later in this section. 

The full polarization function to order 0{a q ) is 



P(q,u) = 



1 



4 y/l - {x + iO) 2 
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The finite imaginary part of P appears when x exceeds 
unity, which we refer to as the absorption threshold. 
Once p\{x) is known, one can use Eqs. @ and (|l7| to get 
the dielectric function e(q,co). For example, in the static 
limit we obtain Eq. ([7]). The real and imaginary parts 
of e(q, lu) as a function of x are illustrated by Figs. [2je) 
and [5Jf) for the case of a suitably low oiq = 0.3. The 
RPA predictions are also plotted in Figs. Me) and [2^f ) 
for comparison. As one can see, the RPA underesti- 
mates Ree(q,oj) at x < 1 and overestimates it at x > 1. 
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FIG. 2. (Color online) Panels (a) and (b) depict function pi(x) where x = uj/v q q. The red dashed line, the blue dotted line, 
and the black solid line are, respectively, the self-energy term, the vertex term, and their sum, which is Pi{x). Although not 
shown in the figure, very close to x = 1 the divergence of the self-energy term overwhelms that of the vertex term, which causes 
the sign change of Kepi(x) at x — 0.998 and Impi(a;) at x = 1.005. Panels (c) and (d) illustrate xpi(x), the quantity which (if 
multiplied by 4i) gives the interaction correction to the universal conductivity, ct/cto — 1. The meaning of the three curves is 
the same as in panels (a) and (b). Panels (e) and (f) depict the real and imaginary parts of the dielectric function for a q — 0.3. 
The green solid line is the RPA result and the black curves (solid and dashed) include the interaction correction. The dashed 
part of the curve is within a q from the absorption threshold. The first-order perturbation theory is expected to be unreliable 
in that region. Unphysical behavior near the threshold is exemplified by the negative sign of Ime in the inset of panel (f). In 
panel (f) the second order interaction correction (black curve) becomes larger than the RPA value (green curve) at u)/v q q ~ 9.1 
(not shown). 



Near the absorption threshold, at |1 — x\ < a q , the first- 
order results are deemed unreliable and are shown by the 
dashed line. In that region higher order corrections be- 
come as important as the first-order ones. An example of 
inapplicability of the first-order perturbation theory near 
x = 1 is the negative sign of Im e(q, uS) at sufficiently large 
a q , see the inset of Fig. |2](f). 

Let us next discuss the analytical results for the func- 
tion P\(q,uj). This function can be written as 



P 1 (g,w) = 2P a (q,Lu) + P b (q,u) 



(18) 



where P a and Pb are the self-energy and vertex terms rep- 
resented by the corresponding diagrams in Fig. [T] (The 
factor of 2 comes from the symmetry of the self-energy 
diagram.) Below we describe these terms separately. 
For the self-energy contribution we have the result 
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where y = u>/(vq) + iO. The expression for function 



I a (y), which is rather cumbersome, is given by Eq. (35) 
(Sec. 151. 



Equation ( 19 ) is written in terms of the "bare" pa- 



rameters. As explained above, we should rewrite it in 
terms of the renormalized ones. To do so we combine 
2P a {q,ui) with the zeroth order polarization function 
P (g,w) [Eq. Q] and get, to the order 0(a g ), 
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The desired renormalized P$(q,oj) and 2P a (q,ui) are, re- 
spectively, the first and the second line of Eq. ( 20 ) . The 



limiting forms of such renormalized 2P a are as follows. 
In the static limit wc find 



2P a ( 9 ,0) = - 
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where G = 0.916 ... is the Catalan constant. In the op- 
tical limit, we obtain 

1 a q 

2P a (q,Lo)~ — ^, uo^v q q. (23) 
16 iuj 

We expect that the accuracy of this expression is im- 
proved if in place of a q one uses ctk, where k is such that 

lj — Vkk. 

Near the absorption threshold, \x — 1| <C 1, we get 



2P a 



J_ Oqq 

16 V a 



21n2- | 

6 



i 



i 



(l-x 2 ) 3 / 2 3(l-x 2 ) 1 /2 



(24) 



Note that the coefficient in front of the dominant (1 — 
x 2-j-3/2 s i n g U l ar ity depends on the renormalization pro- 
cedure. 

Our result for the vertex term is 
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where function h{y) is given by Eqs. ( 37 )— ( 42 ) of Sec. Ill 



and is computed by numerical quadrature. The limiting 
forms of Pb are as follows. In the static limit we find 
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where C b « 0.48. In the optical limit, u> ^> v q q, we have 
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G£« -0.237, (27) 



which is consistent witrPSJ C' b — (8 — 3n)/6. Near the 
absorption threshold \x — 1| <C 1, we find the analytical 
expression 
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which agrees with the result of Ref. [23 further provid- 
ing the subleading divergent term (determined by the 
numerical factor inside the logarithm). 



III. DERIVATION 

In this section we discuss how our results for the self- 
energy and vertex corrections have been derived. Within 
the Matsubara formalism,^ the diagrams we compute 
are expressed by the integrals 
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where ./V = 4 is the spectral degeneracy and the sums are performed over fermionic Matsubara frequencies v = 
-n(2n + 1)0. Green's function G, the self-energy matrix E, and the Coulomb interaction kernel V are given by 



G(k, iv) = (iv -uk-cr) 



(31) 



where <r = {a x ,a y } is the vector of the Pauli matrices. To obtain results at real frequencies uj the analytical 
continuation id — > u + iO has to be done at the end of the calculation, as usual. In the zero temperature limit j3 — > oo, 
the integration over v and the analytic continuation lead to 
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+ £k,q-Ek', q [(l - n k ■ nk+q)(l - ^k' • ftk'+q) + («k x n k+q ) • (n k / x n k / +q )]| . (33) 

These integrals can be simplified by transformation from the Cartesian (k x ,k y ) to the elliptic coordinate system 
(/x, v), where < fi < oo and < v < 2ix. In this system the coordinate grid is made of ellipses and hyperbolas 
with the foci at k = and k = — q. The transformation formulas are k x + ik y = (q/2)[cosh(/.t + iv) — 1] and 
d 2 k = (g/2) 2 (cosh 2 /i — cos^). The integral for 2P a becomes 
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which can be evaluated analytically in terms of the dilogarithm function Li2(z). For oj > we get Eq. (19) with 
2x 2 xh- 2x 2 
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The vertex correction to the polarizability is represented by the integral in Eq. ( 33 1 . Although not amenable to 
analytic evaluation, this integral can still be simplified by employing the elliptic coordinates: 
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x [(sin v sin i/ + sinh/u sinh//) + a; 2 (cosh pcoshp' + tanh p tanh // cos v cos z/)] 



where ?/ = x + iO. Integrating over v and v' ', one is lead to Eq. (251 with If, given by the two-dimensional integral 
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FIG. 3. Integration contours in the complex a plane to eval- 
uate Re/;,(a;) for x > 1 in Eq. (37 1. The dots indicate the 
poles of the integrand. 



the denominator. For such x, one can derive the following 
alternative formulas for Denoting x — cosh(p/2), 
where p > is real, we see that the poles of the integrand 
are at a\,i = p^fb. The result for /(, is complex, with the 
principal value integral (denoted by V) giving the real 
part. The a integral can be computed as the integral 
over contour depicted by the dashed line in Fig. [3j which 
can be deformed into the union of contours C\ and Ci- 
The real part of the integrand vanishes on C2, and thus 
only C\ contributes to Re lb. 

In turn, the imaginary part of If, can be obtained using 
the Sokhotski-Plemelj identity 



where K (z) and E(z) denote the complete elliptic inte- 
grals of the first and the second kind. 

For x < 1, the calculation of Pf, using these formulas is 
easily done numerically and the result is real. For x > 1, 
the standard numerical quadrature routines fail because 
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followed by simple algebraic manipulations. In the end, 
we obtain the alternative formulas for the real and imag- 



the integration path in Eq. (37) passes near the zeros of inary parts of lb suitable for x > 1: 
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(42) 
(43) 
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Straightforward numerical evaluation and asymptotic 
analysis of these expressions lead to the results presented 
in Sec. III1 above. 

Since all the above calculations have been done within 
the Dirac approximation, it is instructive to discuss how 
they can be generalized to a more realistic lattice model. 
We will consider in particular the self-energy term (the 
vertex correction can be analyzed similarly). We will 
show that the correction calculated from the Dirac model, 
Eq. (32), and that computed from the lattice model is 



vanishingly small in the limit of interest, {lu <C v/\a\,q <C 
l/|a|}, where a is the vector connecting a pair of nearest 
lattice sites. 

Let the kinetic energy matrix on the lattice be H = 
"Hk • cr, where %k = (%£,%!<), an d let the on-site den- 
sity distributions be described by a form-factor function 
^(q), which is close to unity at q <C 1/6 and rapidly 
decays to zero at g > 1/6, where 6 ~ a is the charac- 
teristic size of the site orbitals. In terms of these nota- 
tions the self-energy on the lattice E q = E q • <x, with 
E q = (E^, E q ), is given by the equation 



iE q = 



/ 70^(q-k) e ^- k ) °(?^ 



Hi, 



(44) 



Note that the form-factor .F(q) regularizes the short-range behavior of the interaction potential and serves as a cutoff 
on the momentum transfer introduced in Refs. [TBI and IT9l Note also that the integration in Eq. (44) is over the entire 
momentum space to account for Umklapp processes. 

From Eq. (44) we see that the self-energy in the lattice model and the Dirac model [Eq. (11)] have the same 
functional form if the deviation <5k of k from a corner of the BZ K is small: 
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The lattice analog of Eq. ( 32 ) is 
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where the integral is now taken over a Brillouin zone (BZ) 
and N s = 2 accounts for the spin degeneracy only. It is 
easy to see that this integral is convergent and that the 
difference between Eq. ( 32 ) and Eq. 46 is of the order of 
(g/Aiat) ln(g/Aiat). This difference is vanishingly small 
in the limit of interest, as we stated above. 



IV. DISCUSSION 

In this section we compare our work with previous lit- 
erature and discuss some effects beyond the first-order 
perturbation theory. 

The interaction correction to the static dielectric func- 
tion e(q, 0) has been considered in Ref. [30l Our Eq. ^ 
for this quantity differs from the result obtained therein 
by the numerical coefficient of the quadratic term, 0.778 
versus 0.53. This discrepancy is for two reasons. First is 
the apparent error in the numerical evaluation of the ver- 
tex diagram in Ref. 30. Second is the different treatment 
of the self-energy contribution. In Ref. 1301 this contri- 
bution is assumed to be completely absorbed into the 
velocity rcnormalization. In our renormalization scheme, 
only the first term in Eq. ( 19 ) is absorbed while the see- 



not merely a matter of convention. It has to do with the 
principal distinction regarding the relations between ob- 
servable and nonobservable quantities. While the former 
relations do not depend on the renormalization scheme, 
the latter do. In our case, the expansion of the static 
dielectric function e(q, 0) in powers of a q (defined by the 
quasiparticle phase velocity v q ) is unique and given by 
Eq. ([7]). On the other hand, the expansion of e(q, 0) in 
powers of the "bare" coupling a has the form 



e(g,0)~l + |a+ ^0.778 - ^ In ^ 



(47) 



ond leaves a finite remainder, Eq. (21 ). This difference is 



in which the coefficient for a 2 depends on the nonuniver- 
sal cutoff parameter A. Clearly, this coefficient can be 
discussed only after the renormalization scheme is pre- 
cisely defined, as we have done here. 

Next, the optical limit o f the pol arization function has 
been vigorously debatecP^SHiEl] j n con ^ ex t f the 
interaction correction to the universal conductivity of 
graphenc. Our method of calculation and the final re- 
sult, Eq. ([9]), are the same as in Ref. [TBI The interaction 
correction to the optical conductivity of doped graphenc 
has been studied in Ref. 32. In the limit w > /j, where 
\x is the chemical potential measured with respect to the 
Dirac point, one expects to recover the behavior charac- 
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the ladder sum compared to what was obtained by GFM. 
This intriguing problem warrants further study. 
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FIG. 4. (Color online) The real part of the dielectric function 
near the absorption threshold for ot q — 0.3. The black solid 
line is our first-order theory, the red dash-dotted line is the 
ladder sum from Ref. 29 



teristic of neutral graphene. In this limit the results of 
Ref. [32] agree with those of Ref. |T|] and therefore with 
ours. The origin of the discrepancy with Ref. 20] is dis- 
cussed in Appendix |A") 

Finally, the behavior of P(q, to) near the absorption 
threshold \uj/v q q — 1| <C 1 has been studied by the au- 
thors of Ref. [20] (henceforth GFM). In this special re- 
gion the perturbative expansion of the dielectric function 
diverges, see Fig. |4j The divergence of the self-energy 



(1 



-3/2 



[Eq. 



p4[j ] than that 
[Eq. ( 28 ) ] of the vertex term. How- 



term is stronger, 
- (l-a;)- 1 ^^!- 
ever, GFM argued that the divergence of the self-energy 
term can be trivially absorbed into the velocity renor- 
malization thus leaving only the divergence of the vertex. 
Summing the ladder series for the vertex corrections, they 
obtained a nonperturbative expression for e(q, Li). This 
expression is also plotted in Fig. [4] assuming the renor- 
malized velocity of GFM coincides with v q . 

An important qualitative prediction of GFM theory is 
vanishing of Ree at certain x < 1, see Fig. [4j which sig- 
nals the presence of a new collective mode — "excitonic 
plasmon." While we find this prediction very interest- 
ing, we wish to express some reservations in the validity 
of GFM approach. We believe that the nonperturbative 
treatment should begin with the resummation of the self- 
energy not the vertex term because the former is more 
divergent. Such a resummation would make the quasi- 
particlc velocity v q momentum-dependent, cf. Eq. (10 1. 



In other words, the linear Dirac spectrum would be re- 
placed by a spectrum with a finite curvature. As in the 
case of the static response, it is not possible to faithfully 
represent the effect of such a curvature by simply replac- 
ing v with another constant number. It is easy to see that 
the finite curvature of the spectrum modifies the behav- 
ior of the dielectric function over the range of frequencies 
~ aq, which is much wider than the interval a 2 q where 
the higher-order terms considered by GFM are impor- 
tant, see also Fig. [4] Our preliminary analysis suggests 
that this significantly modifies the analytical structure of 



Appendix A: Optical limit 

In this Appendix we review the derivation of the po- 
larization function in the limit v q q -C Li and attempt to 
settle the existing dispute in the literature. In this limit 
the requisite integrals simplify so that they can be evalu- 
ated analytically. Equation ( 32 1 for the self-energy term 



can be expanded to second order in q and reduces to 
f d 2 k 2q 2 (uj 2 +4w 2 fc 2 )E fe 

2P a {q,Uj) ~ J ^ 2 fc 2[ 4u 2 fc 2 _ ( w + l0 )2]2 

The integration is convergent and yields^ 



(Al) 



2P a (q,u) 



1 aq 

16 ioj 



(A2) 



independently of A, which is our Eq. (23 1. In turn, the 



integral in Eq. ( 33 1 for the vertex correction is also con- 
vergent and has the value oP^ 



8 — 3-7T aq 



24 



in agreement with our Eq. (27) 



( |A2| , and |A3]), one arrives at Eq. ^ 



(A3) 

Combining Eqs. 



This straightforward calculation was questioned in 
Ref. [20] (henceforth referred to as JVH) where the twice 
larger value of 2P a was obtained. JVH used the dimen- 
sional regularization (DR) scheme in which the fermions 
are initially assumed to live in D — 2 — e dimensions, 
with the limit e — > taken at the end of the calculation. 
Let us critically analyze this approach. 

Correcting several innocuous mistakes of JVH [such as 
the overall sign error in their Eqs. (62) and its carry over 
into Eqs. (F3)-(F12); also, the missing factor of 4 in front 
of k 2 in the denominator of Eq. (F7)], combining together 
their Eqs. (F7) and (F9), doing the analytic continuation 
from imaginary to real frequency, and finally, temporarily 
setting v = 1, we see that the DR approach predicts 



2P^\q.u) = {D-l) 



d D k 2q 2 {uj 2 +4fc 2 )£D R 



2q z / dkv e {k) 



{2ir) D k 2 [4k 2 - (uj + iO) 2 } 2 
(uj 2 + 4k 2 ) ££ R 



k 2 [4k 2 - (uj + iO) 2 



(A4) 



Here function v e (k) given by 



modified: 



\D-\ 



(47r) D / 2 r(D/2) 



(A5) 



generalizes the free-fermion density of states (DOS) fo(fc) 
to the case of D dimensions and T(z) is the Eulcr 
Gamma-function. The self-energy S}? R in Eq. (A4) is 
given by 



e 

7 



= — A' 



r(i-f)r(^±i)r(^ 

(4ir) D / 2 T(D) 



(A6) 



where we introduced the extra factor A e compared to 
Eq. (65) of JVH, as it is commonly done. This factor 
makes the (energy) units correct and goes away at the 
end of the calculation. However, it is convenient to think 
that A is of the order of uj/v, which is the typical momen- 
tum of electrons and holes that determine the response 
at frequency ui. 



The evaluation of the integral in Eq. ( A4 1 indeed gives 
the twice larger value than Eq. (A2| in the limit e — >• 0. 



Hence, the discrepancy in hand is not due to a mathe- 
matical mistake of JVH. 

The debate held by the authors of Refs. [TS] and UH] 
on the one side and JVH on the other revolved around 
the question of whether the ultraviolet (UV) regular- 
ization is done in a physically correct manner. Exam- 
ples of the regularization schemes that have been dis- 
cussed include the cutoff on the interaction potential, 
V(k) -> V(k) exp(— fc/A) and the cutoff on the band- 
width, Uo(k) — > fo(fc)9(A — k). In the DR scheme, the 
interaction potential V(k) is unchanged but the DOS is 



1 + Q(e) 
X s 



i/ (k) 



(A7) 



Note that both v {k) and v e {k) are the bare DOS not 
the renormalized one, which would follow the spectrum 
described by Eq. ( |To| . 

Th e so ft UV cutoff is assured by the last factor in 
Eq. (A7). It is nearly unity at the physically rele- 
vant momenta k ~ A but gradually decreases with k 
thereby suppressing the bare DOS below that of 2D 
Dirac fermions v$(k). The characteristic momentum cut- 
off scale is A ~ Ae 1 / 5 at which the suppression factor is 
equal to 1/e. The same conclusion follows from the com- 
parison of the self-energy in the DR scheme, Eq. ( A6 ) , 



In 



A 



•0(1) 



(A8) 



with Eq. (111. Accordingly, the comparison between the 



DR and other cutoff schemes can be done by keeping e 
small but finite: 



0(1) > 1 



(A9) 



This kind of interpretation of the parameter 1/e is well 
known in the practice of DR. The problem with the DR 
approach is therefore not the implementation of the UV 
cutoff but the behavior at small k. Within the DR scheme 
the correction to the bare DOS is divergent at small k: 



u (k) 



1 1 A 



(A10) 



Substituting Eqs. (A8)-(A10| into Eq. (A4), we get 



2P Q DR (<z, W ) ~ 2q 2 I dkv (k) (l + ^lnj)^k(c + hiy) - 2U - f ± ^ ~ £ . (All) 

' ' C k J 4 V kj k 2 [Ak 2 - (to + 10) 2 ] 2 8 iw v ' 



In comparison, Eq. ( Al I written in the same units convention reads 

°f p 2 / A\ u? -4- Ak 2 p 2 a 2 





The only difference between the two expressions is in 
the bare DOS. The extra contribution to the integral in 
Eq. (All ) compared to that in Eq. ( A12) comes from the 
product of the artificial v e (k) — vo(Jk) oc \/C correction 
to the bare DOS and the logarithmically large S}? R oc C 



self-energyP2l Since our Eq. (A12) correctly describes the 



physically motivated lattice model (see Sec. Ill I while the 
DR scheme of JHV gives a different result, in our opinion 
the latter does not faithfully represent the behavior of 



electrons in graphene. 

In closing, we would like to comment on the ultra- 
relativistic regime of extremely low frequencies and mo- 
menta at which the renormalized velocity v q approaches 
the speed of light c while a q approaches the fine-structure 
constant ao — e 2 /he sa 1/137. This limit is of no practi- 
cal importance for real graphene samples. According to 
Eq. (13 1, it is reached at unobservably small frequencies 



9 



and momenta of the order of 

A re] = Ae- 4 / Q0 . (A13) 

Nevertheless, it is an interesting limit from a formal 
standpoint because it rep resents a fixed point of the 
renormalization groupP3Sl 

Interaction corrections to the conductivity in the rela- 
tivistic regime have been found to be^ 

f^)_ 1 = ^ ao . (A14) 

(To 187T 

The smallness of the numerical coefficient (92 — 
97r 2 )/(187r) « 0.056, together with the smallness of its 
nonrelativistic counterpart [Eq. ([9])] suggests that such 
coefficient remains small throughout the crossover from 
the nonrelativistic (u) ^> A rc i) to relativistic (uj <C A ro i) 



regime. It is interesting that Ref. [3U also employed a 
DR scheme. Unlike in JVH treatment, the low-energy 
electron DOS has been preserved because the DR was 
performed on the 3 + (1 — e) dimensions of the U(\) pho- 
ton field while the electrons remain in 2 + 1 dimensions. 

Another recent workP^ has presented results for the 
frequency-dependent optical conductivity in the relativis- 
tic regime at the one-loop level of the exact renormaliza- 
tion group on the sa me foo ting as those found in the 
nonrelativistic modeP^ M 20 * 31 * we studied here. We think 
this is misleading for two reasons. First, the two results 
are expected to apply at very different frequency scales, 
u) <C A ro i and uj 3> A re i. Their comparison is meaningful 
only at the common border of validity uj ~ A ro i and only 
by the order of magnitude. Second, it seems that such 
a comparison would require a two-loop calculation, as in 
Ref. M 
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